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Abstract 

The resummation of 0{q^~^^ \t?^ x) terms in the evolution kernels of non-singlet combina- 
tions of structure functions is investigated for both QED and QCD. Numerical results are 
presented for unpolarized and polarized QCD structure functions. 
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1 Introduction 



The resummation of leading small-x contributions to the evolution kernels of QCD singlet struc- 
ture functions [0 may lead to large effects 0. The small-x behaviour of the corresponding anoma- 
lous dimensions is dominated by the leading singularity in the A^-moment plane ~ {as/[N — 1])' 
and corrections to it. For non-singlet kernels such terms are absent both for unpolarized and 
polarized deep-inelastic scattering, and the most singular contributions behave like N^as/N'^Y . 
A resummation of these terms was derived for QCD amplitudes in ref. [§] more than a decade 
ago. 

Similar terms emerge also in QED. There the resummed form of these contributions may be 
described by the structure function method. In the present paper, after setting up our notation 
and recalling the standard NLO formulation in section 2, we present in section 3 a derivation of 
the resummed kernels for the case of QED and QCD |Q in parallel to relate both cases in a direct 
way. The asymptotic kernels are compared with those results found in complete calculations in 
the limit x — *^ up to next-to-leading order (NLO) in QED and QCD. 

Recently very sizeable corrections due to this resummation have been claimed P] for 
the structure functions in both unpolarized and polarized deep inelastic scattering at small 
X. In this way the small-x behaviour of the structure function evolution, as e.g. for 
xF^\x, Q^), f/{x, Q2) _ F2"(x, Q^), and g[{x, Q^) - g^{x, Q^) should be considerably affected. 
In section 4, we perform a detailed numerical analysis and derive the corrections to the various 
QCD non-singlet combination^. 

Also in the case of the singlet anomalous dimensions for polarized deep inelastic scattering 
the leading singularity is expected to behave like ~ N^as/N'^Y . Corresponding equations for 
related QCD amplitudes have been given in recently. The explicit form of the resummed 
anomalous dimension matrix as a function of is derived in where also numerical results 
on the behaviour of the structure function gi{x, Q^) are presented. In the present paper we will 
deal with the different non-singlet cases only and refer for the singlet evolution to ref. 0. 



2 Evolution in fixed— order perturbative QED and QCD 

The evolution equation for the non-singlet combinations q^a,{x,Q'^) of parton densities is given 
by 

The corresponding splitting function combinations P^<^{x,a) are specified below, and ® stands 
for the Mellin convolution, a denotes the running coupling constant in either QED or QCD. In 
order to simplify the notation, we will use the abbreviation a = a/(47r) in the following. The 
scale dependence of the running coupling is defined by 

dn °° 

- E «'^'/^^ • (2) 

fc=0 



rflnQ2 



part of the results has been pubhshed in ref. ji] recently. 
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The first two coefficients of tfie /3-function, Pq and are independent of the renormalization 
scheme. For (one flavour) QED and QCD they read 



QED 



QCD 



4 



^QED 
^QCD 



-4 



34 



C. 



20 



(3) 



CgTrNj - ACfTrN 



f ' 



4/3, Tr = 1/2, and A^/ the number of flavours in 



with Cg = N, = 3,Cf = (N^ - l)/i2N, 
the QCD case. 

In what follows we drop the subscript NS wherever the non-singlet character of the considered 
quantity is obvious from the superscript ±. The splitting functions P^{x,a) are given by the 
combinations 

oo 

P^ix, a) = Pqqix, a) ± Pggix, a) = 5]a'+^P;^(x) . (4) 



1=0 



For the subsequent analysis we restrict ourselves to the consideration of the spacelike case = 
—q^ > 0. Due to fermion number conservation, the expansion coefficients Pf are subject to the 
sum rule ^ 

/ dxPfix) =0, (5) 

JO 

since a acts as an independent parameter. The non-singlet splitting functions for QCD are 
known up to NLO S, P] and read in the MS factorization scheme 



Pqq{x,a) 



2aCB 
+ 



l + x^ 

1 — X 

CIPf{x) + ^-CfCgPg{x) + CfN^TrP^,^.^ 



X] 



—CfCg 
2 



Pa{x) + C(a^ 



(6) 
(7) 



The functions Pi{x), I = F, G, Nf, and A, can be found in refs. 0. The corresponding splitting 
functions for QED in this scheme are obtained from and (|^ by setting Cf = Tr = 1 and 
Cg = 0. Most QED calculations are carried out, however, in the on-mass-shell (OMS) scheme 
for which the NLO splitting functions are different |T^. For x ^ the leading contributions to 
P^{x,a) in the MS factorization scheme are 



p4 


,QED 
-^0 




a) 


= 2a + 2a'^\n^x + 0{a^) 


Px 


,QED 
-^0 




a) 


= 2a - Ga"^ In^ X + 0{a^) , 


p+ 

^ X- 


,QCD 
-^0 




a) 


= 2aCF + 2a^Cl\n^x + 0{a^) 


Px- 


,QCD 
-^0 




a) 


= 2aCF + 2a^ [-3C| + 2CfCg] In^ x + 0{a^ 



(9) 

Since the parton densities q^{x, Q"^) are scheme dependent and hence no observables beyond 
leading order, it is convenient to consider also the evolution equations for the related observables 
directly. These are given by the corresponding structure functions F,^{x,Q'^), obtained by the 
convolution 

Ftix, Q') = cf{x, Q') ® qtix, Q') . (10) 
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Here cf{x,Q'^) denote the respective coefficient functions which can be expanded in a as 

oo 

c,±(x,g2) = 5(l-x) + E«'cJ(x). (11) 



1=1 



After transformation to an equation in a using (H), the evolution equation for F^{x, Q"^) resulting 
from (|l|) and (|l^) reads 



da Pqu"^ 
where in NLO the kernels can be written as 



Kf{x, a) (g) Fj^{x, a) , 



Kfi{x, a) = aPf^sfiix) + 



P±(x)-^PNS,o(a:)-/5o<i(a;) 
Po 



(13) 



The terms oc a(a In^x)*^ emerge in the a-expansion of the kernels Kf(x,a) only in combination 
with the coefficient /So- In this sense the resummation to which we turn now is of leading order. 



3 Resummation of dominant terms in the limit x - 

The most singular contributions to the Mellin transforms of the structure-function 
kernels K^{x,a) at all orders in a can be obtained from the positive and negative 
amplitudes f^{N,a) studied in ref. for QCD via 



M 



K^^oia)] {N) ^ ['dxx'^-'K^_^,{x,a) ^ -]-T^^^{N,a) = ^f^{N,a) 



These amplitudes are subject to the quadratic equations: 



a 


+ 


1 1 r 


N 


87r2 I 


a 
N 


+ 





/o'(iV,«) 



1 1 r 



fo{N,a) 



Here fy{N, a) is obtained as the solution of the Riccati differential equation 



/+(Ar, a) = IQn'ay^ + 2by^^f^{N, a) + 



1 1 



f^iN,a] 



The coefficients and bi in the above relations read for the case of QED 

Oo = 1, ^0 = 1, CLv = 1, by = 0, 

and for QCD || 



do — Cp, 6q — Cp} Civ 



2Nr, 





evolution 
signature 



(14) 

(15) 
(16) 



(17) 



(19) 
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In QED eq. ([17D further simplifies to an algebraic equation with the same coefficients as (0). 
The solutions of (|1^) and (|T6|) were derived in ref. for the QCD caseQ. They are given by 



r;fo''°(iv,a) 



-Nil 



-Nil 



-N ll 



-N I 1 





(20) 



8aCp 



where z = N / ^/2N\^, and Dp{z) denotes the function of the parabolic cylinder 
It is instructive to expand the resummed anomalous dimensions (|20|) and (| 
in /N'^^~^ and transform the result to x-space using 



f211 



into a series 



M 



InM- 

X 



{N) 



k\ 



(22) 



This results in 



K. 



,QED 



X, a) 



2a + 2a^ In^ x + -a^ In^ x + 0{a^ In^ x) 
3 



K^:^Q^{x,a) 
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2a - Qa' In^ x - —a^ In^ x + 0{a* In^ x) , 
3 

2aCF + 2a^Cl In^ x + -a^Cl In^ x + 0{a^ x) 

3 



(23) 



X, a) 



2aCF + 2a^CF 
+ 0(aMn^x) . 



Cf + 



N. 



In^ X + -a^CF 
3 



Cl 



2N. 



In X 



(24) 



Eqs. (p3|-p^ agree with the corresponding result found for Pj 



NS,x^O\-^^ 



x,a) m 



in the complete 

NLO calculations of the non-singlet anomalous dimensions H] in the most singular terms since 



Cn Cf 

^2 



^+'-Cf 



(25) 



holds in SU{Nc). Moreover, K~^^^{x,a) can be compared directly with a result in ref. JTO 
eqs. (2.30, 2.40, 2.43), restricting to the terms oc a^ln^x, where the '-'-non-singlet^ terms were 
given separately in the OMS scheme. Since the corrections there refer to the initial state radiation 
in e~^e~ annihilation the NLO result for a single (massless) fermion line reads 



TV--, QED 
^l,x^O 



OMS 



[x, a) 







2 ^ 






) [ 






2 r 


'2 




) [ 



+ - i In^ X = -6a^ In^ x = KI;^_^q 



MS 



(x,a) . (26) 



^Note that there are a few misprints in eq. (4.7) of ref. 

^Note that the singlet contributions contain terms (x 1/x also in the case of QED. 
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The corresponding result for K^^f^^ can not be derived directly. 



In the evolution equation ([T0| ) aside from the anomalous dimensions -P; (x) also the coefficient 



functions c^l{x] 



contribute. The latter quantities have been calculated to 0{a ) (i.e. 1=2) ||T2| - 
for xF3{x,Q'^) and the non-singlet part of the structure functions F2{x,Q'^) and gi^XjQ"^) 
in the MS scheme. Expanding these coefficient functions for a; — *• one finds after noticing that 
apparent terms oc l/x"^,m = 1,2 cancel in the corresponding expressions of ref. [O, III] 



Q 1 oc as In 



X 



Cj_2 oc In 



X 



(27) 



Hence the terms of O(a^) and O(a^) in ([23H24D can be identified with the parts of the MS 
non-singlet splitting functions oc a(a In^x)'. Thus one obtains 



,QED 



(x, a) 



2,x^0,MS 
p-,QED . X 

p+,QCD_. X 



p-,QCD_, 
2,a;^0,MS 



X, a) 



-a In X 
3 

10 31 4 

a m X , 

3 



10 

y 



-^Cl + ACICg - CfCI ] In^ x . 



(28) 



(29) 



It should be noted that the agreement of the NLO terms between (p^-|2^) obtained from the 
above resummation and (§-|^) holds for < only. This is due to the violation of the Gribov- 
Lipatov relation in the In^ x term of the NLO splitting functions for > 0. 



4 Numerical results for nucleon structure functions 



We now write down the solution of the evolution equation derived in the previous sections and 
then study the quantitative consequences of the leading small-x resummation. In the following we 
will confine ourselves to the QCD case. Here we make use of the fact that the evolution equation 
(p!^ ) for non-singlet structure function combinations reduces to a single ordinary differential 
equation after transformation to Mellin moments. Including the effect of the resummed kernels 
(prp, the corresponding solution can be written as 



F±(iV,a.) = F^{N,ao)(^) 

Van/ 



X < exp 



\ 7NS,o(A^)/2/3o 

da —T^{N, as 



(30) 



2 (3q Jao a 



+ 



as — cio 



2A 







Pi 



itiN) - ^7NS,o(iV) + 2/5oQ,i(iV) 
Po 



with 



and oo 



-2 / dxx 
Jo 



N-l 



[X 



cHN) 



(31) 



In (|30|) 7^ (A^) stands for the two-loop anomalous dimension 'jtiN) with the 



leading 1/A^^ term obvious from (|^) removed^. This latter contribution is already included in 



In (|30| ) and ( p2D we have corrected two trivial misprints in cqs. (21) and (23) of ref. 
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the exponential factor, which in turn is connected to (|2l|) via the subtraction of the contribution 
hnear in a^, 

r±(iV, a,) = r±_^o(iV, a,) - ^ hm [iV7NS,o(iV)] = T^^oi^, a,) + . (32) 

The well-known NLO evolution of F^{N,as) is entailed in (|30|) by simply expanding the expo- 
nential to first order in and Qq. Finally, the transformation of the solution back to x-space at 
any x and affords only one standard numerical integral in the complex A^-plane |15 . 



The remaining quadrature in (p(J| ) can be performed analytically for the '-|-'-case, resulting in 

Jao 2 Oo VOs Cto/ i 0,s ttQ J 

_NA (1- Vl-Aa,)(l + yi-Aao) , ^ 8CV 

2 ''(1 + Vl-Aa,)(l- Vl-Aao) ' iV2 ■ ^ ^ 

On the other hand, the corresponding integration has been carried out numerically for the 
combinations involving the parabolic cylinder function Dp{z = N / yflW^g). Alternatively, one 
can expand the resummed kernels T^{N,as) and r~(A^, a^) in the strong coupling a^, using the 
Taylor series of the square root and the asymptotic expansion[] of Dp{z), respectively. In the 
practical applications considered below, one finds that, even at the lowest x- values considered, 
more than 90% of the resummation effects in (pOD arise from the first two terms beyond NLO in 
the expansion of r^(iV, a^). 



Let us consider the quantitative consequences of the resummation (|21| ) for two representative 



non-singlet combinations. For this purpose, we choose Qq = 4 GeV^ as our reference scale 
in (PPD, and employ the same initial distributions F^(A^, oq) and Aqcd for the NLO and the 
resummed calculations. The evolution is performed for Nf = 4 active (massless) quark flavours. 
Unless another value is stated explicitly, we take A = Ajj^{Nf = 4) = 230 MeV in 



1 



(3i lnln(QVA^ ^ 
'M ln(QVA2) 



1 - 



(34) 



We start with the unpolarized case, where we investigate the evolution of the '+'-combination0 

F2''{x, Ql) - F/"(x, Ql) = c%{x, Ql) ® \ [xu, - xd, - 2{xd- xu)] (x, Ql) , (35) 



adopting the input densities from the MRS (A) |JT6[ global fit. The small-x behaviour of the 
most relevant quantities is given by xu^{x,Ql) ~ x^^^^, xdv{x,Ql) ~ x°-'^^. Note that these 
distributions are rather 'steep', i.e. their rightmost singularity in the complex A^-plane lies about 
0.5 units or more to the right of the leading singularity of the non-singlet splitting functions 
at A^ = 0. Studying the evolution of this F2 difference at very small x is mainly of theoretical 
interest, since it is orders of magnitude smaller than F2^ and F^"' there. In figure 1 the result 
of this investigation is depicted down to x as low as 10~^^. Even at these extremely low values 
of X, the effect of the resummed anomalous dimensions stays at the level of 1% or below, and is 
still dominated entirely by the first two terms beyond NLO. 



''For X <10 ^ the asymptotic series is no longer reliable and one has to refer to the resummed result directly. 
*For corresponding results on the '-'-combination 0:^3 (x, Q^) the reader is referred to ref. 
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In the polarized case we consider the corresponding difference 



9i'{x, Ql) - gn^, Ql) = c-,{x, Ql) ® - (Au, - A4) (x, Ql) . (36) 

This case is practically much more interesting, firstly since - unlike in the unpolarized case - 
the non-singlet distributions are not a priori suppressed here with respect to the singlet ones 
at very low x, and secondly since the shapes of the polarized initial distributions are not yet 
well established. We illustrate the strong dependence of the resummation effects on the latter 
quantities by choosing two partly rather different input sets for Au^, and Ad^. At first, we 



take those of CW |jT^, using xq = 0.75 in eq. (12) of ref. [|T^], which have been used in several 
theoretical investigations 0, The small behaviour of this input is relatively fiat, Au^ ~ 



and Adv ~ ^+0.29 small x. As an example for a more recent parametrization we adopt 
the 'standard' NLO set of GRSV as an input, using their value of the scale parameter. 



^Ms(^/ = 4) = 200 MeV. Here one has Au^ ~ Ad^ ~ x~^'^''', hence the steepness is 

similar to that of the unpolarized initial distributions above. We have evolved both distribution 
sets in NLO from their respective input scales, = 10 GeV^ in [17| and = 0.34 GeV^ in 



l8| , to our reference scale Ql = 4GeV^. 



Before we derive the quantitative results, we notice that eq. (30) violates the fermion number 
conservation for the '-' non-singlet combinations. Here the conjecture is that the coefficient 
functions cfi{x) do not contain terms oc In^ x in the MS scheme. For this no proof exists 
yet, however, we have verified this behaviour up to 2-loop order in section 3 for the coefficient 
functions of xF^, F^^ , and gi^. One should recall that the main resummation effect comes 
from that and the next order. Under this assumption fermion number conservation has to be 
restored for r~_^Q(iV, a^). We approach this problem in several ways numerically. In a first set 
of calculations we subtract a corresponding term oc 6{1 — x) from the kernels K~ derived from 
(|2T|) , in each order in a^. In A^-space this prescription (denoted by 'A' below) leads to 

r-(iV,a,) ^ T-{N,as)-T-{l,a,) . (37) 

Another possibility is the restoration of fermion number conservation by subleading 1/A^ pole 
terms. An especially simple choice (denoted by 'B' in the following) is to modify according 
to 

T-{N,as)^T-{N,as)-{l-N). (38) 

Besides these two prescriptions, which are analogous to the procedure in the second reference 
in [0 with respect to energy-momentum conservation in the unpolarized singlet case, we will 
also show the results for two other assumptions, namely ('C') 

T-{N, a,) T-{N, a,)-{l-2N + N^) (39) 

and ('D') 

r-(iV, a,) ^ T-{N, a,) ■ {1 - 2N + N^) . (40) 

Clearly, the results of the resummed calculation are only trustworthy, and this approach is to 
be preferred over a fixed order calculation, if the difference of the results obtained by all these 
procedures is small. 



The corresponding results are presented in figure 2. For the relatively fiat CW input |T7| , 
the effect is up to about 15% at x = 10~^. However, in the kinematical range accessible for 
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polarized electron and proton scattering at HERA |jT9|] it again again amounts to about 1% or 
less. Note, moreover, the wide spread of the results in dependence of the employed fermion- 
number conservation prescription. Obviously the resummed contributions do not sufficiently 
dominate with respect to subleading terms at any foreseeable energy. For the steep GRSV input 
the effect is of approximately the same marginal size as in the unpolarized case considered 



above. In figure 2 relative corrections are shown. Recall that the absolute values for gi obtained 
in the different parametrizations extrapolating from the range x > 10~^ of the current data down 
to smaller x values vary strongly M,M. 



5 Conclusions 

We have investigated the resummation of terms of order a'^^ In^' x, derived in ref. for the 
QCD case, on the small- a; behaviour of non-singlet functions in QED and QCD. The comparison 
with the corresponding contributions obtained in the same order by complete NLO calculations 
shows the equivalence of both approaches in this limit up to order in both QED and QCD. 
Since the coefficient functions up to two-loop order for the non-singlet combinations considered 
contain only terms less singular in In x in the MS scheme, the contributions oc In^ x in the 
three-loop MS splitting functions P^{x, a) have been predicted on the basis of this resummation. 

A numerical analysis has been performed for the QCD case of deep-inelastic (polarized) 
lepton scattering both off unpolarized and polarized targets. It turns out that the all-order 
resummation of the terms 0{a''^^ In x) leads only to corrections on the level of 1% in the 
unpolarized case of F2 — F2 even down to extremely small x values, x = 10~^^. The corrections 
can be larger in the polarized case, up to about 15 % at x ~ 10~^, depending on the presently 
not yet well established small-x behaviour of the polarized parton densities. In any case, the 
resummation effects are on the level of 1% in the kinematical range accessible experimentally at 
present or in the foreseeable future. 

Presently unknown terms which are suppressed by powers of Inx in the splitting functions do 
contribute in a potentially significant way to the evolution even at the lowest x- values considered. 
This has been demonstrated for the '-' combination gf — g^ by applying several prescriptions 
to implement fermion number conservation into the resummed evolution equations. Moreover, 
the resummation corrections are dominated by the first two terms in an expansion beyond 
NLO. All this indicates that fixed-order perturbation theory remains the appropriate theoretical 
framework for the evolution of non-singlet structure functions even at very small x. 
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Figure 1: The small-x Q^-evolution of the unpolarized non-singlet structure function combination 
— in NLO and the absolute corrections to these results due to the resummed kernel 
derived from ref. 0. The initial distributions at Ql = A GeV^ have been adopted from |jl6 . 
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Figure 2: The relative corrections to the NLO small-x Q^-evolution of the polarized non-singlet 
structure-function difference Qi^ — Qi^ due to the resummed kernel. 'A', 'B', 'C, and 'D' denote 
the different prescriptions for implementing the fermion number conservation discussed in the 
text. The initial distribution at Qg have been taken from ref. [0 in (a) and ref. in (b). 
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